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Abstract. There are few data available with which to constrain the thermal 
history of the intergalactic medium (IGM) following global recombination. 
Thus far, most constraints flow from analyses of the Cosmic Microwave 
Background and optical spectroscopy along a few lines of sight. However, 
direct study of the IGM in emission or absorption against the CMB via 
the IS hyperfine transition of Hydrogen would enable broad characterization 
thermal history and source populations. New generations of radio arrays are 
in development to measure this line signature. Bright foreground emission 
and the complexity of instrument calibration models are significant hurdles. 
How to optimize these is uncertain, resulting in a diversity in approaches. 
We discuss recent limits on line brightness, array efforts including the new 
Large Aperture Experiment to Detect the Dark Ages (LED A), and the next 
generation Hydrogen Reionization Array (HERA) concept. 



1. Introduction 

Characterization of the intergalactic medium (IGM) and early generations of 
luminous compact objects before global reionization (z > 6) is a frontier in 
observational cosmology. For the first ~ 100 Myr following recombination and 
decoupling of radiation and matter, evolution was dominated by gravity acting 
on diffuse distributions of dark and baryonic matter. Models are robust owing 
to linearity of interaction. For the remainder of the first ~ 1 Gyr, nonlinear 
processes become increasingly relevant, e.g., with the widespread star formation, 
and understanding is limited. 

There are no data with which to directly constrain evolution of the IGM 
prior to formation of luminous objects (a.k.a. the "dark age," despite diffuse 
radiation). Even during reionization of small pockets initially and large expanses 
later, there are few data available to constrain the thermal history of the IGM, 
source populations, and how these formed. An integrated global estimate of the 
redshift of reionization assuming a sudden transition (z = 11.0 ± 1.4) stems 
from measured large scale polarization of the cosmic microwave background 
(CMB, [1]). Recent observations of the kinetic Sunyaev-Zeldovich effect on the 
CMB anisotropies have given a model dependent constrain on the duration of 
reionization Az < 7.9 with reionization beginning at Zfj eg < 13.1 at 95% CL ([2]). 

Limits on the neutral Hydrogen (HI) column have been inferred from optical 
and near-infrared (NIR) spectra of luminous objects, where a trough blueward 
of Lya appears for z > 6.3 (the "Gunn- Peterson Trough"). Neutral fraction is 
difficult to infer over a wide range of z owing to the high absorption cross section 
and ready saturation. However, inference may be drawn from modeling of the 
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Lya damping wing, including the effective radius of the proximate ionized zone. 
As of this writing, only two quasars are known for z > 6.4 ([4]), the most distant 
at z ~ 7.09 ([3]). Troughs have also been observed toward gamma ray bursts, 
and there are 3 known for z > 6.4, up to z ~ 8.3 ([5], [6] and [7]). Modeling of 
spectroscopic data for quasars and GRBs enables inference for a handful of lines 
of sight. Broad sky coverage may be limited by nature - there are relatively few 
supermassive black holes that evolved in the first ~ 1 Gyr - and instrumental 
reasons - wide-field NIR surveys and high sensitivity spectroscopic follow up is 
difficult (and subject to atmospheric windows.) 




Figure 1. Average sky brightness temperature of the A21cm line (i.e., spin 
temperature x optical depth) up to the end of reionization [9] . Variation with redshift is 
driven by varying coupling between radiation (the CMB) and gas kinetic temperatures. 
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Figure 2. Power 
in fluctuations as a 
function of redshift 
for four spatial fre- 
quencies on the sky 
(adapted from [10]). 
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The A21 cm transition of Hydrogen at high redshift has not yet been detected, 
but as it would trace the diffuse IGM horizon to horizon without the need for 
illuminating compact objects, study of this transition has the potential to provide 
a new avenue by which to infer the evolution of the IGM and source populations. 
This tracer complements CMB studies and optical/NIR spectroscopy and 
spectral-line origin in principle allows discrimination by redshift, possibly into the 
Dark Age. The physics of the A21 cm transition has been described in reviews 
(e.g., [8]; [10], and references therein), but we summarize here salient points, 
referring to Figures 1 and 2. 

Following recombination, the IGM cooled (adiabatically) faster than the 
CMB. When there was effective collisional coupling between A21cm transition 
and gas temperature, the line will be seen in absorption against the CMB (both 
in the sense of a sky average and fluctuations). This occurred for z S> 40, 
when gas density was high, but before the first stars. Coupling faded with 
expansion of the Universe, reducing contrast between the HI line and CMB. 
With early star formation, leakage of Lyman photons into the IGM recoupled 
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Table 1. High Redshift A21cm Experiments 
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^ Shaded denotes efforts to measure the angle-averaged HI spectrum. LEDA and 
DAWN share the LWAl aperture. DAWN applies beam forming on cold patches of the 
sky. LEDA instead cross-correlates the aperture (§3.5). ' 2 - ) IN: Pune, India; NL: Exloo, 
Netherlands; WA: Boolardy, Western Australia; ZA: Karoo, South Africa; US: VLA 
site, New Mexico. ' 3 - ) Approximate collecting area, mid-band, per polarization. ' 4 - ) 
Maximum array baseline. Redshifts omit a 10% guard band. 

the spin and kinetic temperatures via the Wouthysen-Field effect, resulting in 
a second absorption signature (e.g., [9]). The absorption depth and redshift of 
the minimum depend on the extent and evolution of heating sources (e.g., star 
formation, X-rays from black holes, shocks), which may involve "exotic" physics 
(e.g., dark matter annihilation). Eventually, IGM heating drives the HI signature 
into emission before the HI component is largely destroyed; in this regime, the 
sky-averaged and fluctuations signal will be more extended in fractional redshift 
and weaker in magnitude. 

2. Observations of the A21 cm line 

Sensitivity to the anticipated unpolarized mK signal will be limited by 
systematics. These will arise principally from errors in instrument calibration 
and foreground models subtracted from the data. The two will be derived from 
similar data and likely to be coupled. Unpolarized foregrounds will be smooth- 
spectrum mix of extragalactic point sources and Galactic diffuse emission with 
angular scales overlapping those anticipated for the A21 cm signal ([12]). 

Foreground brightness well exceeds the A21 cm signal. The former rises to 
O(10 3 ) K at 100 MHz even away from the Galactic plane (proportional to 
approximately v~ 2 ^ [13]). This is 0(1O 4-5 ) times the A21cm signal in sky- 
averaged power at ~ 70 and 140 MHz (Figure 1). Fluctuations in foreground 
emission for "quiet" patches of sky are 0(1-10) K) ([11]) and may be only 
O(10 2 " 4 ) times the A21 cm fluctuations of (9(1-10) mK) at ~ 140 MHz (Figure 2). 
The brightness of individual foreground sources extends down to the A21 cm 
signal. This can be subtracted where it is above limits imposed by source 
confusion and contamination by sidelobes. However, this exceeds thermal noise 
limits for even modest integrations with arrays matched to the large angular 
scales of fluctuations. Below this nonthermal limit, other approaches are required 
(e.g., subtraction of an extant high-resolution sky or a statistical model). The 
limit reported by [11] is a few mJy at 150 MHz and ~ 2' resolution (~10K). 

A more detailed study of foreground structure is still necessary for EoR 
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observations. At present, prediction of the unpolarized diffuse component 
requires interpolation between widely spaced frequencies. Diffuse emission 
dominates the polarized sky (cf. pulsars and compact AGN). There are no wide- 
angle low-frequency polarized sky surveys or consistent picture from the few 
patches that have been studied, and angular and frequency structure in Stokes 
Q and U can differ considerably from that of Stokes I ([12]). 

The point source population is somewhat better understood. Catalogs of 
low-frequency point sources with arcminute resolution are available from the 
Cambridge surveys at 151-178 MHz (e.g., [16]) and from the VLA 74 MHz survey 
([17]). In the Southern Hemisphere, Culgoora Surveys at 80 and 160 MHz ([18]) 
provide measurements at arcmin resolution. A GMRT survey that will provide 
finer resolution at 150 MHz 1 is underway. 

For single dipole instruments, a careful instrumental design coupled with a 
polynomial fit has been demonstrated to be effective in foreground subtraction 
([19]). However, uncertainty in antenna gain patterns (e.g., in response to 
the environment) and unmodeled responses of electronics have thus far limited 
sensitivity to the global 21cm signal other than for sudden reionization ([19]). 

For dipole arrays, a higher degree of calibration is possible, benefitting efforts 
to detect A21 cm fluctuations, or total power via array elements capable of 
measuring calibrated total power. Bright compact sources are used to calibrate 
direction and frequency dependent instrument response, and they may be 
subtracted directly from visibility data ([41], [50], [11], [22], [23]). Foregrounds 
are mapped in parallel and, given a smooth and slow variation with frequency, 
galactic diffuse emission and sub-confusion level point sources can in principle 
be filtered ([24], [25], [26], [27], [28]) provided that calibration errors do not 
result in pixel to pixel spectral fluctuations (e.g., [29]). Equivalently, calibration 
errors can couple a fraction of the polarized emission from Galactic sources into 
the unpolarized visibilities (as noted above), contaminating the A21cm signal 
([11], [30]). For idealized models of sky and instruments, it has been possible to 
filter out the polarization leakage by identifying the effects of different Faraday 
screens along various lines of sight and subsequently filtering the data ([31]). 
Active areas of investigation are in correction of ionospheric distortions ([32]), 
wide field polarization ( [33] ) , correction for antenna gain patterns ( [34] , [35] ) and 
direction dependent deconvolution ([36], [37]). 

3. Ongoing A21 cm Experiments 

Arrays thus far reflect a diversity of architectures, which is a consequence of 
incomplete knowledge of the low-frequency sky, complex coupling of dipole 
elements to the sky (and man-made transmissions), and uncertainty regarding 
optimization of foreground subtraction. The RF hierarchy (i.e., layout of 
antennas) dictates the challenges to be met in signal processing and in calibration 
and imaging. Equivalently, the available computing budget where signal 
processing is implemented with general purpose platforms (e.g., BlueGene, 
GPUs) constrains what RF hierarchies are practical. We briefly review 
ongoing experiments (Table 1), emphasizing differences in architecture and signal 
processing. Recent limits on the A21cm signal are presented in Figure 3. 

1 http://tgss.ncra.tifr.res.in 
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Figure 3. Experimental constraints on reionization thus far. (left)- Lower limits 
for the breadth of reionization (Az) for a "sudden reionization" model, estimated from 
singlc-dipolc total-power measurements by EDGES ([19]). Shading denotes 68% and 
95% C.L. (right)- Upper limits on 3D power spectra (z=8-9) from Wcsterbork data (*) 
([11]), and the GMRT for two different foreground subtraction techniques (adapted from 
[38]). Linear fits over 2.0(D) and 0.5 MHz (o) provide a foreground subtraction model. 
The 0.5 MHz limit may preclude a cold IGM at z~9 (dot-dashed line; [26]). 



3.1. Giant Metrewave Radio Telescope (GMRT)-EOR experiment 

Fourteen GMRT dishes are located in a 1 km 2 area, providing the low surface 
brightness sensitivity required to attempt detection of the A21 cm line. A 16 MHz 
wide 150 MHz receiver band corresponds to z~8.5. The narrow ~3.3° field of view 
(FOV; half-power) of the dishes has enabled the experiment to begin without 
wide-field polarization calibration and correction for ionospheric distortions. 
The target field is centered on a pulsar and in a novel calibration scheme, 
gated visibilities (pulsar ON - OFF) are differenced, thus eliminating the other 
(unswitched) foregrounds and simplifying phase and polarization calibration 
([39]). Foreground subtraction has thus far depended on fitting a frequency linear 
baseline on 0.5 and 2 MHz scales in visibility space, leveraging the anticipated 
smooth spectrum of the foregrounds ([38]). 

3.2. Low Frequency Array (LOFAR)-EoR experiment 

Twenty-three high frequency (120-240 MHz) LOFAR stations are arranged in 
the ~3km core for use by the EOR experiment. Each station comprises two 
clusters of 24 tiles of 16 dual polarization dipoles (768 total). This is the deepest 
hierarchy of RF-element clustering for any EOR instrument. Two levels of beam 
forming within the LOFAR stations enables good signal rejection but at the 
cost of complex sidelobe structure that may be difficult to map in assessment 
of coupling to the sky brightness distribution (cf. the GMRT). Calibration 
and visibility plane foreground subtraction will be executed via generalized self- 
calibration techniques that incorporate direction dependent effects ([40], [41]). 
This is supported by stations far from the core, which are used to catalog point 
sources and to accumulate enough lines of sight to enable 3D tomography of the 
ionosphere ([42]). The EOR experiment is in its first year of data collection. 
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3.3. Murchison Wide-field array (MWA) 

The MWA will comprise 128 tiles of 16 dual polarization dipoles, each with 
a beam former establishing a single ~20° half-power FOV at 150 MHz. EOR 
observing will utilize a strongly centrally condensed pseudo-random configuration 
of ~ 1 km 2 ([43]). Dense u,i>-plane sampling and co-planarity will enable imaging 
via co-adding of warped snapshots ([33]). As such, wide-field correction can be 
achieved with geometric precision, and ionospheric distortion (approximated as 
a 2D rubber sheet) can be implemented in image space. The large number 
of baselines will enable detection of many calibrators per snapshot, a running 
estimation of the many, variable tile gain patterns, and correction during gridding 
of visibility data. Source subtraction down to nearly the confusion limit will be 
accomplished via peeling in visibility space and image-based Forward Modeling 
will enable deconvolution ([36], [37]). Among EOR projects, so great reliance on 
image-based processing is unusual, but if successful, it may be relevant to much 
larger arrays, for which the volume of of visibilities, ocO(N^ ipole ), may preclude 
traditional techniques. 

3.4- Precison Array to Probe the Epoch of Reionization (PAPER) 
PAPER comprises 64 broadband dipoles (100-200 MHz) with plans to build out 
to 128 or 256 ([44]). A shaped ground screen for each reduces what would 
otherwise be horizon-to- horizon dipole response to ~60°. LOFAR, MWA, and 
PAPER represent a progression to flatter hierarchies vis-a-vis RF and signal 
processing architecture. PAPER operates in a drift scan mode and benefits from 
the stability of the dipole gain patterns (cf. MWA and LOFAR). The array can 
be reconfigured, enabling experimentation with pseudo random distributions for 
building a sky model and redundant distributions that can benefit calibration 
and concentrate sensitivity on particular angular scales ([45]) as may be desired 
for first detection of A21 cm brightness fluctuations. 

3.5. Large- Aperture Experiment to Detect the Dark Ages (LED A) 
LED A 2 is distinct among arrays in that it will attempt detection of the A21 cm 
signal at z ~ 20 (30-88 MHz band) at the end of the Dark Ages and thus to 
constrain initial conditions for reionization (e.g., [9]). The target is the angle- 
averaged signal, and a large-N array is used enable joint estimation of instrument 
calibration and a sky model (in contrast to what is possible with a single dipole). 
LEDA will attempt to leverage the O(10 x) magnitude of the A21 cm absorption 
trough at z~20 compared to emission at lower redshift. A 512-input full-Stokes, 
FPGA/GPU-enabled correlator (e.g., [46]) will be integrated into the lOOxllO-m 
Long Wavelength Array (station 1) for 30-88 MHz observing. Outrigger dipoles 
will be instrumented for total-power measurement using calibrated noise sources. 
Their gain patterns will be calibrated using cross-correlation data for baselines to 
the core stations. The outriggers will be well separated from the core to reduce 
mutual coupling effects that contribute high-order terms in gain patterns, to 
provide baselines that resolve most diffuse Galactic emission, and to improve 
sensitivity to the point source population. 



2 http://www.ledatelescope.org 
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4. Hydrogen Epoch of Reionization Array (HERA) 

The Hydrogen Epoch of Reionization Array (HERA) is a planned 2 nd generation 
(mid-decade) experiment that targets the detailed evolution of the A21 cm power 
spectrum (e.g., [47]), and to a lesser extent imaging ionized "bubbles" around 
selected luminous early quasars 3 . In concept, HERA will have similar sensitivity 
to phase 1 of the Square Kilometer Array, though HERA will be an experiment 
optimized for study of the A21 cm line at z > 6. It will be timed to follow current 
generation US projects (PAPER, LEDA, MWA) and exploit lessons learned. 

The HERA conceptual design comprises a centrally condensed, filled aperture 
with a sparsely sampled "halo" of outlying antennas. In seeking to characterize 
A21 cm power spectra at a detailed level, emphasis is on constraint of systematics 
during instrument design and in foreground characterization and removal. The 
primary specifications affecting sensitivity are FOV, pass band, total collecting 
area, and area per element (Table 1). The geographic distribution of elements 
enters vis-a-vis the quality of foreground mapping and source subtraction. FOV 
and pass band dictate directly the range of /c-space to which the array is 
sensitive. For large k, i.e., fine angular scales, artifacts from subtraction of 
point-sources may become increasingly serious ([49]). Errors in reconstruction 
and subtraction of diffuse foreground emission, particularly in the presence of 
ionospheric fluctuations, will also reduce sensitivity over a range of k, although 
quantitative modeling of the effect has been limited. 

Figure 4. Indicative 3-D power spec- 
trum measurement error vs a model, for 
a filled-aperture 10 5 m 2 array, assuming 
a neutral IGM, lOOO' 1 integration, 8MHz 
bandwidth, T sky = (470, 1250) K, and 
5000 antennas in a filled core and 3km 
diameter envelope. Symbols indicate in- 
dependent bins. The vertical dotted line 
corresponds to bandwidth in k units. The 
effects of calibration error and systemat- 
ics are ignored. Adapted from ([48]). 

The LOFAR architecture presents one model by which to accumulate 
collecting area (a sparse array with beam forming of many-dipole stations). In 
contrast, HERA is anticipated to use a "flatter" array configuration. This can 
be more computationally challenging, starting with cross-correlation, which is an 
0(N 2 ) problem. For a collecting area of 10 m 2 per antenna per polarization (more 
than for PAPER and LAVA, and half that for MWA), 10 4 elements are required, 
for which correlation of 100 MHz demands 40 PFlops -1 . A scalable solution 
for correlation drawing on GPU computing has been proposed by [46], with 
demonstration planned by LEDA. Calibration carries a similar computational 
order (e.g., [50], [51]), and scalable demonstration relevant to HERA may be 
achieved by LEDA and or MWA. 

Summary— We have reviewed efforts to detect the high redshift A21 cm 
signal with radio arrays. Diversity of approach reflected in the design of first 

3 http: / /reionization. org/RFI2_HERA.pdf. This described a roadmap of instruments. Here we 
refer to the next, mid-decade experiment, dubbed HERA-2, simply as HERA. 
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generation devices marks an "epoch of experimentation" from which strategies 
and optimizations for next generation arrays, such as HERA, will be based. 
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